Few-photons model of the optical emission of semiconductor quantum dots 
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The Jaynes-Cummings model provides a well established theoretical framework for single elec- 
tron two level systems in a radiation field. Similar exactly solvable models for semiconductor light 
emitters such as quantum dots dominated by many particle interactions are not known. We access 
these systems by a generalized cluster expansion, the photon-probability-cluster-expansion: a reli- 
able approach for few photon dynamics in many body electron systems. As a first application, we 
discuss vacuum Rabi fiopping and show that their amplitude determines the number of electrons in 
the quantum dot. 
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Introduction: For describing the quantum dynamics of 
semiconductor light emitters involving weak light-matter 
coupling and large photon numbers, an expansion involv- 
ing mean field quantities and their fluctuations (often 
called cluster or correlation expansion) provides a well 
controlled theoretical scheme to treat many particle sys- 
tems of electrons, phonons and photons How- 
ever, for semiconductor emitters operated in the limit of 
few photons 5] and few electronic levels, such as quantum 
dot- wetting layer (QD-WL) systems (cf. Fig. [l]i)) [1], 
the cluster expansion (as illustrated below) breaks down, 
since it is based on the assumption that fluctuations of 
mean field quantities have minor influence. With the rise 
of quantum information and the search for proper solid 
state emitters, QD-WL systems for single and entangled 
photons are in the focus of current research [1, 0, 0| ■ 
Therefore, there is an urgent need for the extension of 
the standard cluster expansion. Such an extension - the 
photon probability cluster expansion (PPCE) is devel- 
oped in this Letter. 

We generalize the cluster expansion to the limit 

of few photons in systems with few quantum confined 
electronic levels such as semiconductor QD-WL systems, 
cf. Fig. [T^). The strength of our method is also to keep 
the accurate results of the standard cluster expansion for 
large photon numbers 0, 01 and the strong coupling 



a) 



HzzT 

WL QD WL 

J=L 




Figure 1: Scheme of the model system: a) a QD is embedded 
inside a WL leading to confined QD states coupled to a WL 
continuum, b) the emission of the confined states couples to 
a photon mode in a cavity. 



limit for high quality cavities 

To illustrate the standard cluster expansion and its 
break-down, we introduce c^, c as the Bosonic construc- 
tion/destruction operators for a single photon mode with 



frequency ujq and 



as the Fermionic operators for 



electrons. Subscripts c, v, denote the lowest confined 
conduction and valence band state of the QD, respec- 
tively. The corresponding model of the QD-WL sys- 
tem is depicted in Fig. [T] b). In the following, this 
model will be used to illustrate our approach to repro- 
duce exactly the Jaynes-Cummings model (JCM) (one 
electron approximation in atomic two level systems (i)) 
and extent it to many electron systems (QD-WL sys- 
tems (ii)). To start, we illustrate the difference between 
the single and many electrons systems: The Jaynes- 
Cummings model (JCM) Q, containing only one elec- 
tron in an atomic two level system is exactly solvable by 
expanding the solution into states having distinct pho- 
ton number and electronic state. Therefore, quantities 
of interest such as photon number (c^c), photon-photon 
correlation (c'c'cc), electronic densities {ala^), (ala^) 
can be calculated exactly. Here, the expectation value 
is denoted by (...) — tr(. ..p), where p is the statis- 
tical operator. In particular, by the exact expression 
{ala^alac) — (ajoc) for a single electron, the hierar- 
chy of quantum correlations is closed Q. Unfortunately, 
however the JCM does not provide a scheme to include 
corrections resulting from further fermionic particles be- 
yond the two level, one electron approximation. Beyond 
the JCM, i.e. in many electron systems, such as an oc- 
cupied WL, Fig. [1] non-factorizing correlations such as 
(ala„ala^)'' — (alayalac) — {a\ac) ^ are important 
[j, |lO| . This is the case in particular for QDs embedded 
in a WL [4] , since in this case the electrons in the WL and 
in the QD are Coulomb coupled: Electrons are scattered 
from the WL to QDs and vice versa, which lead to a cor- 
relation between electrons in the QDs and in WlQ]. The 
first correlation contribution requires an electron and a 
hole to be combined for emission. 
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So far the standard cluster expansion of (ala^ala^)'^ 
provides access only to large photon number dynamics or 
in weak coupling-short time limit {ti\g'"'^\-t <C 1 , where g"'^ 
is the light matter couphng constant, see below). This 
is, because the standard cluster expansion constitutes a 
hierarchy of equations for correlations, in which the dy- 
namics of relevant observables are not connected directly 
to the photon number state. Examples include quanti- 
ties such as {ala^c^c), which include contribution from 
photon states with one and more photons . This will be 
illustrated later using Eq. , it is clear that this kind of 
expansion is not suitable, if correlations with a distinct 
number of photons are of crucial importance. 

A theoretical scheme, which attacks this problem for 
a single photon mode, but is easy to extend to multiple 
photon modes, is presented in this Letter. This new ap- 
proach will allow to treat few photon dynamics for quan- 
tum information in many particle systems such as QD- 
WL systems with few electronic levels, that contribute to 
the emission. 

Model system: We illustrate the proposed photon- 
probability cluster expansion (PPCE) using the basic 
light- matter Hamiltonian H = Hq + Hei^pt'^- 

Hq = hwQC^c + hecCila^ -f he^ala^ (1) 
Hei-pt = -hig'^'^alaj +g''-*ala,c), (2) 

(illustrated in FigUlD)). For simplicity, we consider only 
a valence band v and conduction band c state with the 
energies fe„ and fee respectively. The off-diagonal cou- 
pling matrix g^'^ denotes the electron-photon-interaction 
strength. Additionally, the two level QD is assumed 
to be embedded, in an electronically occupied WL. In 
such a many particle system the JCM-valid factoriza- 
tion {ala^ala^) = (aja^) has to be replaced by 
(ala^ala^) = {ala^){l - (ajaj) + {ala^ala^)" [Ij 
Clearly, because the one electron assumption in the QD- 
WL system fails, Pauli-blocking terms are introduced [j]. 
Since the correlation (a'^a^ala^)'^ cannot be treated ex- 
actly, an analytical solution of the model is not possible 
anymore, and the cluster expansion of {ala^a\a.^) pro- 
vides only a perturbational approach. 

We start showing that the cluster expansion yields un- 
physical behaviour in the few photon strong coupling 
limit for two electronic levels: Up to the order Ig"'^^, 
using the Heisenberg equation of motion, the cluster ex- 
pansion reads 3, 4]: 



5t(ctc> 



= i{ev - £c + Wo + i-ip){ala^c^) - ig""" ftp 
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(4) 



Here, we introduced the hole occupation — 1- 
the electron occupation = {a^a^ and the phenomeno- 
logical pure dephasing constant 7^. The spontaneous 
emission term /J^ is equal to in case of the JCM and 
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Figure 2: The JCM (solid/red), standard cluster-expansion 
atom like using one electron assumption (dashed/green), and 
standard cluster-expansion using /Jp = f'^f'^ in QD-WL sys- 
tem (dotted/blue) and PPCE using f!p^„ = in QD-WL 
system in new scheme (dotted dashed/black). With the ini- 
tial conditions: = 0.95, po = 0.8, pi = 0.2 and a 
coupling constant \g^'^\ = 0.5 1/fs for a QD resonant with the 
cavity mode. ( (a)7p — l/(10ps) and (b)7p — 0/(ps)) 



f^p = (ala^ala^) = ff" for the QD-WL case 0. In 
Fig. [2] we plot the dynamics of the electron occupation 
for: (i) exact JCM model [8], (ii) cluster-expansion using 
one electron assumption /|p = f^, (iii) for the QD-WL 
system /|p = f^f^ and also (iv) for the below developed 
PPCE, introduced later. 

The exact JCM solution for the one electron assump- 
tion leads to well known vacuum Rabi flops. Fig. [2]shows 
clearly that the standard cluster expansion of the JCM 
and the QD-WL system lead to negative and unphysical 
results in the few photon limit. However, we note that 
in the many QD and/or many photon case (c^c) 3> 1, 
correlation expansion provides an excellent tool 0, [lo| . 
As expected, all approaches agree for small times (weak 
coupling limit h~^\g"'^\t ^1). 

The breakdown of the cluster expansion results from 
the fact that the system is completely characterized with 
mean single particle quantities. This is problematic in 
the limit of few electrons and photons, since correlations 
(in particular fluctuations) between electrons and photon 
number states are crucial in this limit. To improve this 
situation, we focus now on the development of the PPCE. 

Photon- Probability Cluster Expansion (PPCE): Impor- 
tant quantities to characterize a quantum optical field 
are the normalized intensity-intensity correlation func- 
tion g^'^\T = 0) = ^'^etcp^ ^'^'^ ^^'^ photon intensity 
g(i)(r = 0) = (c'l'c) [9|. Both quantities can be measured 
and they are directly connected to the photon number 



and intensity-intensity expectation value 



The idea - to circumvent the problems characteristic for 
the cluster expansion (cf. Fig. [5]) - is to formulate all 
quantities in terms of n photon probability p„ = ( |ri) (yij) , 
where \n) denotes the Fock-state of n photons in the sys- 
tem. This idea is already succesfully applied in JCM [81], 
but needs to be generalized to the many electron case. 
The PPCE introduced now combines the advantages of 
the exact JCM-solution and the standard cluster expan- 
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sion. The observables read in the Fock basis: 

OO CSC 

(c^c) = ^npn, (c^c^cc) ^^n{n- l)pn. (5) 

Eq. ([5]) show that it is sufficient to calculate the equa- 
tions of motion for the probability to find n photons 
Pn = in the system. Heisenberg equations for 

Pn, using the Hamiltonian, Eqs. dHl]), couple to assisted 
occupation probabilities: (|n) (njata^), (|n) (rt|aj,a^) and 
transitions (|?i)(n-f l|aja^) or similar quantities if a more 
general system is considered. We discuss the application 
of the theory to the JCM (one electron) and the QD-WL- 
system (many electrons) separately: 

For the case of one electron {JCM case), the equation 
of motion for the probability p„ reads: 

dtPn - -2V^lm{g--{\n){n-l\ala^)) 

+2\Amim(5''^(|n-f l)(n|aK))- (6) 

The n-photon probability Pn is driven by the photon 
transition assisted polarization (|n -I- l){n\al^a^y. 

dt {\n + I) {n\ala^) 

= ^(£^, - £c + t^o + ^lp){\n + l){n\ala^) 
+*5-*V^(p„+i - ./^i) - i9^'^*V^rn.{7) 

Besides the free energy rotation quantum optical absorp- 
tion and emission processes occur in the last line. These 
processes are related to /j^ = Pn — (|n)(n|aja„) and 
fn = (l'^)('^kc'^c) for hole and electron densities assisted 
by n photons. Finally, the electron and hole densities 
fn'^ couple to quantum optical transitions: 

cJj;^/'' = 2V^Im(g-(|n+l)(n|ata,)). (8) 

To close the hierarchy of equations of motion we used 
(a|a^afca/|7i)(n|^ = strictly valid only for the single 
electron case, since two electrons are annihilated. Eqs.(IB]- 
[5]) without pure dephasing reproduces the JCM Q . This 
is an important benchmark of the new PPCE. 

If we consider more than one electron [QD-WL case), 
dt{\n + l){n\a\ja^ , couples to (ajaja^a„|ri)(n|) ^ 0. 
In order to truncate the PPCE hierarchy of equations 
of motion at least approximately, we need factorization 
rules for correlation between electrons, but also for other 
bosons to e.g. introduce cavity loss and radiative decay 
consistently. 

We use and generalize two basic factorization limits, 
Hartree-Fock and Born approximation (i,ii), respec- 
tively. 

(i) The Hartree-Fock approximation (ala'jakai) « 

(a|a;) (ajofc) — (a|afc)(a]ai) is based on the as- 
sumption [l|, [l^, that the density matrix p can 
be described by single particle quantities using the 
Generalized Canonical Statistical Operator (GCSO) 



p{t) « exp[-^^^. A,j(t)aJa^./(fcBT)]/Z ^ with the La- 
grange parameters (t) . For the n-photon number con- 
tribution of the density matrix {n\p(t)\n) /pn, a photon 
number dependence of the Lagrange parameters is intro- 
duced: {n\p{t)\n)/p„ « eM-E;j Kjit)aUj/ikBT)]/Z. 
From this, we recognize the factorization rule: 
{\n) {n\a\ala3a4) = {{\n) {n\a\a4) {\n) (nlalas) - 

{\n) {n\a\a3){\n){n\ala4))/pn. 

(a) The second order-Born approximation for a corre- 
lated electron-photon system embedded in a bosonic bath 
(e.g. phonons or other photon modes , b) assumes that 
the density matrix can be factorized: p pb ® Pei-pt, 
which implies: (^\n) {n\a\a2b\b4) ~ (^\n) {n\a\a2) {blb4) . In 
particular, for QD-WL system, to some electronic states 
bath like properties (e.g. WL states) can be attributed, 
so a factorizing density matrix p = Pei-uncorr Pph+corr 
can be assumed. This leads to a factorization rule: 
(|n)(n|a-aj) w {\n){n\'){alaj) for WL states i,j. 

Our factorization rules can be extended to include 
high order PPCE correlations, preserving the important 
electron-photon number correlations. We now evaluate 
the equations of motion in the PPCE approach: In 
comparison to the single electron case, the equation for 
Pn (Eq. ([5])) does not change. Eq. ([7]) is modified for the 
QD-WL system case using the modified Hartree-Fock 
factorization rule: 

dt {\n + l){n\ala^) 

= i{sy - ec + uJo-\- «7p)(|n + l){n\ala^) 

+«5''Vn+ 1(P«+1 - fn+i){Pn+l - fn + l)/Pn+l 

-*.g-V^/,^/,yp„. (9) 

In comparison to the one electron assumption, Eq. ([7]) 
the terms {pn-i-i — fn+i) and are replaced with {pn+i — 
fn+i)(Pn+i- fl+i)/Pn+i and fl'J^/pn respectively. Sim- 
ilarly, Eq. dH]) does not change. (In Eq. ([9]), some 
sourceless terms including {\n){n\aXac) are neglected.) 
Eqs. dH),®,® constitute the new quantum optical dy- 
namics for a few photon/ few electronic level system. 

Discussion: The best benchmark for the new PPCE 
scheme is the case of vacuum Rabi flops where quan- 
tum fluctuations dominate the dynamics. For strong 
coupling, the in-scattering rate from the WL has to be 
insigniflcant compared to g^^. Ref. If suggest for g'"'^ 
chosen here, WL densities <C lO^^cm"^. Now, we reit- 
erate our discussion of Fig. [2l Obviously, only the PPCE 
(dotted dashed) provides physical reasonable results as 
opposed to the standard cluster expansion method, since 
it guarantees positivity of the electronic density. In par- 
ticular, we recover the exact JCM solution. The Pauli 
blocking term occuring in the QD-WL system does not 
lead to strong deviations for not long times from the JCM 
case. For the calculations in Fig. [21 we assumed, that 
initially there is one electron {f^ -t- (1 — /'*)) = 1 in the 
QD = f^. For the numerically evaluation, we choose 
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Figure 3: The PPCE solution for semiconductor QDs. With 
the initial conditions: = 0.3 and = 0.1, po = 0.05, 
pi = 0.95 and a coupling constant g — 0.5 meV for a QD 
resonant with the cavity mode and 7p — (10ps)~^. In (a) 
the different lines are (solid/red), (green/long dashed), 
g^(0) (blue/short dashed), (c) The maximum value of (?'^' 
for PPCE is plotted over the number of electrons in the QD, 
with pi — 1.0 (black/straight) and p2 — 1.0 (red/dashed). 

as initial correlations a finite photon number distribution 
Pn(io)) whereas the initial transitions (|rt+l) (n|aja^) will 
be set to zero, i.e. assuming an uncorrelated initial state. 
Also vanishing initial correlations between electron and 
hole system are assumed: The photon number, electron 
and hole densities will be set to = • Pn{to) and 
fn = •Pn(*o)- Comparing the PPCE calculation to 
the cluster expansion and the JCM, only the PPCE re- 
covers the exact JCM result. Therefore, we assume that 
it is more valid for the QD-WL case than the standard 
cluster expansion. The damping of Rabi flopping de- 
scribes dephasing due to electron-phonon coupling and 
Coulomb correlations betweeen QD and WL |1L] . 

We now focus on the specific QD-WL properties. In 
Fig. [3l the dynamics of /'', g'-^-' and 5^^^ is plotted for 
semiconductor specific initial values ^ , occuring in 
QD-WL systems The average photon number oscil- 
lates between and 1, the g^^^ between and 0.7 is in 
the anti bunching regime 0] . After the first Rabi flop the 
system goes into the single photon emitter regime with 
a g^'^^ < 0.5[T3|. In contrast to the JCM, the strength 
of the vacuum Rabi flopping is reduced in dependence on 
how strong the electron and hole density deviates from 
the one electron assumption = f^. This implies that 
the amplitude of the Rabi flops might be used as a mea- 
sure for the number of electrons in the actual QD. Since 
the amplitude of the Rabi flopping is directly connected 
to g^^\ it may be possible to stcar g'-^'> with the devia- 
tion from the one electron assumption. This is illustrated 
more clearly in Fig. [3t), where the number of electrons 
in the QD is plotted as a function of the maximum of 
(7^^^ , which is directly connected to the maximum of vac- 
uum Rabi flopping. It is seen that a maximum of g^^^ 
is achieved, if we have only one electron fe = fh- This 
behaviour can be attributed to the Pauli-blocking terms 
fe ' fhJ because they drive the emission and are enhanced 
in the one electron case. 



In conclusion: a modified cluster expansion scheme 
to describe, semiconductor quantum dot - wetting layer 
quantum optical devices in the strong coupling limit in- 
volving few photons and few dominant electronic states 
is introduced. This photon-probability-cluster expansion 
(PPCE), reproduces the Jaynes-Cummings model (JCM) 
as a well known benchmark in quantum electrodynamics, 
the standard cluster expansion cannot. For the QD-WL 
system, we have shown that the amplitude of vacuum 
Rabi flopping and the maximum of g^^^ depend on the 
numbers of electron and holes pumped into the quantum 
dot device. This is a first important prediction from the 
PPCE. 

We envision, that the new scheme can be applied to 
many solid state structures with few electronic states 
limit in the few photon limit. The approach will stimu- 
late further investigation beyond JCM: Our approach al- 
lows to include electron-phonon interaction and Coulomb 
interaction additionally to the already demonstrated 
electron-photon interaction in a similar fashion. 
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